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In the last years a general consensus has emerged on the use of ultra-high energy cosmic rays 
(UHECR) data as a powerful probe of the validity of special relativity. This applies in particular 
to the propagation of cosmic rays from their sources to Earth, responsible for energy suppressions 
due to pion photoproduction by UHE protons (the Greisen-Zatsepin-Kuzmin limit) and photodis¬ 
integration of UHE nuclei (the Gerasimova-Rozental limit). A suppression in the flux of UHECRs 
at energies above 40 EeV - as expected from both these interactions - has been established ex¬ 
perimentally beyond any doubt by current experiments. However, such an observation is still not 
conclusive on the origin of the suppression. In particular, data from the Pierre Auger Observatory 
can be interpreted in a scenario in which the suppression is due to the maximum acceleration 
energy at the sources rather than to interactions in the background radiation. In this scenario, 
UHECR data can no longer yield bounds on Lorentz invariance violations which increase the 
thresholds for interactions of nuclei on background photons, in particular through modification 
of the dispersion relations. Here we argue in turn that the study of UHECRs still represents an 
opportunity to test Lorentz invariance, by discussing the possibility of deriving limits on violation 
parameters from UHECR phenomena other than propagation. In particular we study the modifi¬ 
cations of the shower development in the atmosphere due to the possible inhibition of the decay 
of unstable particles, especially neutral pions. 
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1. Introduction 

Our universe appears to be very well described by quantum mechanics at small scales and 
by general relativity at large scales, but the correct way to completely unify these two theories is 
still unknown. It is possible that the correct theory of Quantum Gravity (QG) would predict that 
space-time is subject to quantum fluctuations and the geometry of the world, with all its symme¬ 
tries, emerges in the semiclassical limit. Lorentz Invariance (LI) is therefore not guaranteed to be 
an exact symmetry, and violations (LIV) can therefore be possible. Although these violations may 
only be very small, since no measurement has as yet found any evidence, it has been shown in the 
last two decades that measurable deviations can be present even at energies much lower than the 
QG scale. In particular, possible LIV effects could show themselves in Ultra High Energy Cosmic 
Ray (UHECR) phenomena. 

Although the idea of possible experimental effects of EIV is quite old, the possibility of putting 
extremely strong limits on, at least some, EIV parameters from UHECRs detection was firstly 
quantitatively discussed in [1] and later on refined in many ways. It has to be recalled that UHE¬ 
CRs are the most energetic particles in the Universe, and the use of the experimental observations to 
constrain possible EIV effects is natural, in particular since it has been shown [2, 3] , (the Greisen, 
Zatsepin and Kuzmin - GZK - effect) that these particles do not propagate freely in the Universe 
but interact with the photon backgrounds (Cosmic Microwave Background Radiation, CMB, and 
ExtraGalactic Background Eight, EGBE) that fill it. The processes involved, photopion production 
and photo-dissociation, implicate low energies at terrestrial accelerators. On the contrary, they ap¬ 
pear to occur at UHE given the small energy of the background photons. 

Consequently, as soon as the evidence of the suppression in the spectrum of UHECRs around 
5 X 10^^ eV became undisputable, based on results from HiRes [4], Auger [5] and, more recently. 
Telescope Array [6], limits on those violating parameters were derived. A discussion and refer¬ 
ences can be found in [7]. All these bounds however rely on the assumption that the observed 
spectrum is generated in sources which inject nuclei at energies larger than those characteristic of 
these interactions. 

Here we discuss the status of these bounds in the light of the combined analysis [8] of the measure¬ 
ments of the Pierre Auger Observatory both of the spectrum [5, 9] and nuclear composition [10] 
which strongly suggests that the observed suppression in the spectrum might be due to the maxi¬ 
mum cosmic ray acceleration energy at the sources rather than to an effect of their propagation in 
extragalactic space. This fact was recently noted by several authors [11, 12, 13, 14, 15]. 


2. LIV and UHECR propagation 

In recent times possible El violating models have been deeply analyzed and compared with 
available experimental data [7, 16]. In particular, in the approach of Effective Eield Theories, pos¬ 
sible Eorentz (and CPT) violating terms have been described in the Standard Model Extension 
(SME) in [17]. When symmetry (and renormalizability) is no longer a guide, the number of pos¬ 
sible terms is essentially infinite. Those that can be subject to experimental verification (several 
hundreds) are described in [16]. 
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To parametrize departures from relativistic invariance, as for instance in [1], we only consider 
modification of the dispersion relations (which in the SME approach corresponds to modification 
of kinetic terms in the Lagrangian), which amounts to assuming that the relation, connecting the 
energy and momentum of a particle, is modified as: 

Ef -pj= ^ + ^£2+^ ^2.1) 

where p = \~^\, p is a. function of momenta and energy, Mp ss 1.2 x 10^^ eV, the Planck mass, 
is the possible scale where QG effects become important and ft, which can be either positive or 
negative, parametrizes the strength of LIV for particle i. The last equality reflects the fact that LI is 
an exceedingly good approximation of the physics we know, so that modifications are expected to 
be quite small, making an expansion of the LIV dispersion relation in terms of E/Mp appropriate. 
Here we make the simplifying assumption that all possible terms at the same order have the same 
form, which amounts to neglect higher order corrections to Eq. 2.1. Also, we consider only the 
lowest correction term in E/Mp present (higher orders are totally negligible at UH energies). In 
practical terms, only n= 1,2 will be relevant [1]. 

The right hand side of Eq. 2.1 is invariant when fi = 0. We will assume normal conservation of 
energy and momentum. In a LIV theory it is always possible to move all the violations to the 
dispersion relations [18]. Einally we assume that, in nuclei, LIV only affects nucleons: since then 
the relevant energy entering the modified dispersion relations is the energy/nucleon, this implies 
that violations are milder for nuclei than for protons; for a nucleus of atomic number A, effectively 
Mp AMp. 

Lrom Eq. 2.1 it is clear that the correction term is always much smaller than both {E^,p^) even for 
E « 10^® eV. However, as soon as' p > the correction becomes larger than the 

mass of the particle, and this can lead to very important effects [1]. We consider here how LIV 
affects the thr'eshold energy for the GZK process pYbkg —^ {p^n)n, where Y},kg is a photon of the 
CMB or EGBL. The threshold for this process, in a LIV world, is modified: 


Egzk 


mpfUn 


Egzk 


p{Ep,pp,mp,Mp)p{E,^,Pjc,mn,Mp) 

2C0y 


( 2 . 2 ) 


(cOy being the energy of the background photon). The last equation has to be solved for Ep = Egzk- 
Lor our simplified treatement, we will assume that /, are the same for all the hadrons. 

The case we are interested in is for /,• < 0. As soon as /, moves from zero towards negative 
values the threshold energy at first slightly increases, but for fi < -2.5 X 10-'^ {n = 1) [/; < 
—4 X 10^^ {n = 2)], Eq. 2.2 has no longer real solutions [1]: the photo-pion production reaction is 
no longer kinematically allowed and protons propagate freely in the Universe, if only interactions 
on the CMB are taken into account. Interaction on the EGBL are affected for larger values of —fi?' 
Lor nuclei, for which the relevant process of interaction on the universal backgrounds is photo¬ 
disintegration, an equation corresponding to Eq. 2.2, with Mp —)• AMp, can be written. The modifi¬ 
cation of the thresholds is similar to that for protons. 

'since at the leading LIV order E ~ p we will use them without distinction. 

^For positive fi the thresholds move lo lower values. However other processes are allowed in this case, such as 
Vacuum Cherenkov p —> py, for which very strong bounds exist [19], from the mere existence of UHECRs. 
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Figure 1: Simulated energy spectrum of UHECRs (multiplied by E^) at the top of the Earth’s atmosphere 
with maximum source rigidity Rmax = 5 x 10*^ V and 7 = 2, along with Auger data points. The propagation 
is simulated switching off the interactions with photon backgrounds. Partial spectra are grouped according 
to the mass number as follows: A = 1 (red), 2 < A < 4 (grey), 5 < A < 14 (green), total (brown). Right: 
average (top) and standard deviation (bottom) of the X^^ax distribution (assuming EPOS-LHC for UHECR-air 
interactions). Gray bands show the systematic uncertainties. 


Limits on LIV parameters derived from the observed steepening of the spectrum of UHECRs have 
been reported in literature [20, 21]. These limits, however, depend crucially on the assumption 
that the observed flux suppression is originated by the propagation of UHECRs. Auger compo¬ 
sition data combined with those on the all-particle spectrum have been analyzed in [8] (see also 
[11, 12, 13, 14, 15]) and they indicate, at least when interpreted in a simple source model, a differ¬ 
ent scenario in which the flux suppression reflects the end of cosmic ray acceleration at the source, 
with the implication of requiring hard injection spectra. 

It is therefore worthwhile to verify if ElV can still be bound in this scenario. 

To simulate ElVs we have propagated UHECRs switching off the interactions with background 
photons, only accounting for energy losses due to the expansion of the Universe. To account for 
these losses we used both a simplified version of SimProp [22] and a semi-analytical code, within 
a simple source model consisting of identical sources, uniformily distributed in comoving volume 
throughout the whole Universe, emitting all nuclei in a rigidity dependent way. These nuclei are 
the same impinging on the atmosphere with energy reduced by the redshift. Eor this analysis we 
assume injection of 4 nuclei (H, He, N, Ee). The absence of interactions might make source evolu¬ 
tion effects important. We have checked that this does not happen for soft injection spectra as the 
one discussed here. 

It has been shown [22, 13] that it is not possible to reproduce the data with a single extragalactic 
component, a different one being needed below E = 5 x 10^^ eV, unless possible primary interac¬ 
tions are considered in the source environment [14, 15]. Eor our analysis we will only consider 
larger energies. 
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We arbitrarily impose fixed 7 = 2, consistent with Fermi acceleration mechanism, and maximum 
rigidity Rent = 5 x 10^^ V, and search the best composition that reproduces the experimental data, 
taking into account only statistical errors. The results for spectrum and composition are presented 
in Fig. 1: the resulting composition is light to intermediate (H, He, N). The overall (spectrum and 
composition) reduced « 2.3 is acceptable given the simplifications made and does not allow us 
to exclude LIV in the propagation of nuclei. 

It is however obvious that the above statement cannot be taken as evidence of LIV, since many 
other astrophysical/particle physics explanations can be considered. 


3. LIV effects on shower development in the atmosphere 

In principle, all aspects of UHECR physics can be modified by LIV. For instance, LIV can 
affect the cosmic ray acceleration processes, and also the energy losses during acceleration. As an 
example, consider LI violation with /, > 0. In this case, as soon as the energy of the accelerated 
nucleus is above some threshold, vacuum Cherenkov process becomes possible and essentially no 
further acceleration is possible. This case has been discussed, in a different context, in [19]. 

On the other hand important effects are expected in the interactions of UHE particles in the atmo¬ 
sphere and in the decay of secondary particles. These effects can make some parts of the kinemat- 
ical space forbidden for the processes and therefore make some reaction impossible. 

Consider the most important decay for atmospheric showering, —)■ 77 . It can be shown that, 

considering EIV modihed dispersion relations, the kinematics of the decay changes into ([23]): 

'4 + ^(A4+”-/7(4+"+4^"))-2(Er.^o.-FriF:.) = 2FriFn(l-cos0i,2) (3.1) 

Since there are very strong limits [24] on fy we will assume it to be zero. The right hand side 
of Eq. 3.1 is non negative, while the left hand one can become negative for large enough Ejc and 

<0. Therefore neutral pions do not decay if x 10*^ eV {n= 1) 

(« L4/|/;t|*^‘^ X 10 '^ eV, n = 2). The case with n = 2 is expected to give a very small effect and so 
it will not be considered in the following. 

To test this effect we have generated 10^ atmospheric showers with a modified version of CONEX 
[25] in which the decay of unstable particles is inhibited for E,- > (Mpm^/I/,-1)1/3 fixing = 1 . ^ 

The results of this simulation are presented in Eigs. 2. In the left panel, the expectation for X^ax vs 
energy for El shower development (solid lines) and EIV case (dashed lines) is reported. As one can 
see, the change in the particles’dec ay (mainly neutral pions) has the net effect to move the shower 
maximum to higher altitudes as the electromagnetic part of the shower consumes faster. Erom the 
observational point of view the EIV effects on the shower development make the interactions of 
nuclei (and protons) primaries compatible with (El) interactions of heavier nuclei. 

We incidentally note that another effect of the suppression of neutral pion decays at high energies 
would be the increase of the number of their interactions. This would give more muons at ground 
level with respect to El shower development, as shown in Eig. 2 (right panel) where the average 

^For the case considered, the effect of LIV can be opposite for particles with respect to antiparticles. This effect has 
been neglected in the present work. 
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Figure 2: Left panel: Expectation for {X^ax) vs energy for LI shower development (solid lines) and LIV 
case (dashed lines). Right panel: Average number of muons vs primary energy in LI and LIV cases. 



Figure 3: Same as in Fig. 1. In this cas a maximum source rigidity Rmax = 10*®'^ V and y = 2.65 are used. 
Here EPOS for UHECR-air interactions has been modified so as to take into account the LIV effects. 


number of muons vs primary energy in LI and LIV cases is reported. A study of possible constraints 
to LIV from observational data on muons in showers will be addressed in a forthcoming analysis. 
The modifications of the shower development due to LIV (Fig. 2, left panel) do affect the analysis 
of the spectrum and composition data. Including such modification, the solution reported in Fig. 1 
corresponds to a reduced 2 ^^ > 4. However a better description of the data is found at 7 = 2.65, 
Rcui ~ 1.6 X 10*^ V with a reduced ~ 1-6. In this case, reported in Fig. 3, the composition 
changes and is dominated by light elements (H, He): there is in fact no need of heavier primaries, 
given the change of the interactions in the atmosphere. It has to be reminded, however, that the 
modihcations of the shower developement described here are only relevant for the case of rather 
strong violations i.e. n = 1 and /,■ « 0(1). Milder violations give smaller effects. Similar effects in 
a LI context were presented in [26]. 
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4. Conclusions 

UHECRs are the highest energy particles since almost immediately after the Big Bang. It is 
therefore natural to use their propagation and interactions to probe the structure of space-time, and 
in particular possible tiny violations of Lorentz Invariance. Experimental data, both on the spec¬ 
trum and chemical composition, from the Pierre Auger Observatory have been recently analyzed 
and the suppression of the all particle spectrum has been interpreted as a possible manifestation of 
the maximum energy reached at sources that accelerate nuclei in a rigidity dependent way. 

We have analyzed the status of limits on LIV in this scenario, assuming violation in the prop¬ 
agation through the universal background radiation: for values of the violation parameter /, < 
—2.5 X 10-14 {fi < —4 X 10 ’), interactions on photons are switched off, still experimental data 
can be reasonably described, so that essentially no bounds on LIV parameters can be cast. 
Violations can also affect the interaction of UHECRs in the atmosphere and in particular the shower 
development. We have found that even in this case we cannot exclude this effect, that is however 
much weaker than that relevant for propagation. The results reported here cannot be taken as evi¬ 
dence for LIV. Even for the hard injection spectra astrophysical explanations are available (see for 
instance [27] and references therein), and (non LIV) changes in the interactions at the highest ener¬ 
gies are possible. However, finer spectrum and composition measurements at the highest energies 
can conceivably be used to put constraints on (at least some) violations of relativistic invariance. 
Einally we remark that we have discussed here only a small number of cases of all the possible vio¬ 
lations. These are probably the most important affecting UHECR propagation and interaction, but 
certainly a general analysis of all aspects of LIV on cosmic ray physics would be very welcome. 
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